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Abstract 

Graphene, which as a new carbon material shows great potential for a range of applications because of its exceptional electronic 
and mechanical properties, becomes a matter of attention in these years. The use of graphene in nanoscale devices plays an impor- 
tant role in achieving more accurate and faster devices. Although there are lots of experimental studies in this area, there is a lack of 
analytical models. Quantum capacitance as one of the important properties of field effect transistors (FETs) is in our focus. The 
quantum capacitance of electrolyte-gated transistors (EGFETs) along with a relevant equivalent circuit is suggested in terms of 
Fermi velocity, carrier density, and fundamental physical quantities. The analytical model is compared with the experimental data 
and the mean absolute percentage error (MAPE) is calculated to be 11.82. In order to decrease the error, a new function of E 
composed of a and P parameters is suggested. In another attempt, the ant colony optimization (ACO) algorithm is implemented for 
optimization and development of an analytical model to obtain a more accurate capacitance model. To further confirm this 
viewpoint, based on the given results, the accuracy of the optimized model is more than 97% which is in an acceptable range of 
accuracy. 
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Introduction 

The astonishing discovery of graphene as an extraordinary two- 
dimensional (2D) material with low dimensional physics, and 
possible applications in electronics [1-6] has attracted the atten- 
tion of scientists in these days. Geim, in 2004, demonstrated 
that the six-membered rings are the basis of all carbon ma- 
terials in electrochemical biosensor research [7]. The remark- 
able electrical properties of graphene such as fast electron trans- 
port, tunable energy bandgap, high thermal conductivity, and 
ballistic transport at room temperature give rise to the potential 
applicability in electrolyte-gated transistors [8-11]. Graphene, 
as a nearly perfect 2D crystal free of the structural defects 
[12,13] shows ballistic transport because of its significant high 
electron mobility at low temperatures, which can reach up to 
200,000 cm 2 /V-s with a typical carrier concentration of 
2-10 11 cm -2 [7,14]. Recently attempts have also been made to 
use graphene as a novel channel material in field effect transis- 
tors (FETs) for electronics [15]. The remarkable properties of 
graphene reported so far included high stiffness with a Young's 
modulus of approximately 1 000 GPa, a significant heat conduc- 
tivity of 3000 W-(m-K) -1 , and large specific surface area of 
2600 m 2 -g _1 [15-17]. Intrinsic graphene is a semi metal or a 
zero band gap semiconductor, which results in a high electron 
mobility at room temperature [18]. Therefore, the electron 
transfer in graphene is expected to be 100 times faster than that 
in silicon. Other advantages of graphene, which make it a 
perfect semiconductor is its massless Dirac fermion structure 
with zero band gap (graphene is considered to be theoretically 
lossless) [19]. Compared to silicon-based devices, graphene 
with its outstanding properties such as consuming less energy 
and faster heat dissipating show a great promise in electrolyte- 
gated graphene field-effect transistors (EGFETs) [20]. 

An EGFET fabricated on a SiCVSi substrate with gold source 
and drain electrodes and a graphene layer as a conducting 
channel can be seen in Figure 1. A 300 nm SiC>2 layer as a 



back-gate dielectric has been deposited above the doped silicon 
substrate. The graphene layer, the gate, and a quasi-reference 
electrode were covered by a small droplet of ionic liquid [21]. 
The standard three-electrode electrochemical cell using a poten- 
tiostat has measured the interfacial capacitance of the graphene 
[22]. Vq is the voltage applied at the Pt gate electrode, and V re f 
is the voltage measured on the quasi-reference electrode. Ye et 
al. have discussed the distinguished properties and high 
carrier-density transport of ion-gated mono-, bi-, and trilayer 
graphene bases in double-layer transistors [23]. They demon- 
strated that V re f ~ 0 over the whole sweep range of Vq, which 
has led to dropping all the applied Vq at the liquid/graphene 
interface [15]. 

In order to minimize the background capacitance,the mentioned 
configuration is employed which can also prevent the graphene 
edges from exposure to the electrolyte. To interpret the elec- 
trical response of the device, an analytical model of the EGFET 
together with the equivalent circuit describing its operation is 
discussed in this paper. As depicted in Figure 2, the measured 
capacitance is assumed as the contribution of two interfacial 




Electrolyte-Gated Field Effect Transistors Equivalent Circuit 

Figure 2: A cross-section of graphene-based electrolyte-gated field 
effect transistor, together with the equivalent electrical circuit. 
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capacitances which arise from the double layer formed by ions 
at the graphene-ionic liquid interface and the quantum capaci- 
tance of graphene. The particular case of EGFET is discussed in 
the context of 2D systems. 

To substantiate all the consideration just made, the equivalent 
circuit of the EGFET device is suggested in terms of a simple, 
well-defined theoretical model. It is seen that Cq is the capaci- 
tance that forms between the gate and the ionic liquid, is the 
electrical resistance of the solution, Cq represents the geomet- 
rical capacitance of the double layer/graphene interface to 
model the accumulation of a layer of counter-ions on a charged 
electrode. Finally, Cq is the quantum capacitance of the EGFET 
associated with the finite density of states of graphene [24]. 
Figure 2 shows that Vq has a strong influence on the capaci- 
tance. The total capacitance is given by 1/C = I/Cq + 1/Cq with 
the smaller of the two capacitances dominating the total capaci- 
tance. Previous experimental studies have reported a large 
geometrical capacitance (several tens of uF/cm 2 ) [25,26]. Since 
the two mentioned capacitances are connected in series, the 
smallest one would dominate the total capacitance. Hence, 
geometrical capacitance is neglected compared to the theoretic- 
ally predicted quantum capacitance of graphene. As expected, 
Cq dominates the total capacitance, which is why the position 
of the Fermi energy E? can be tuned by applying only small 
values of Vq. These explanations of the current study are 
consistent with those of Jilin Xia [26] in 2009, who found that 
the Debye ionic screening length of the ionic liquid is virtually 
zero, which makes the quantum capacitance a dominant source 
of the measured capacitance. They performed the measurement 
of quantum capacitance of bilayer graphene in an ionic liquid 
electrolyte. The aim of this study is to evaluate the quantum 
capacitance of single layer graphene sheet as a function of 
voltage, and validate theoretical predictions with the experi- 
mental results [26]. 

Results and Discussion 
Proposed model 

The quantum capacitance of nanoscale devices is considered as 
an important quantity in the design of nanoelectronic devices. 
The classic expression for quantum capacitance utilized in the 
prediction of the theoretical model for an ideal single layer 
graphene [27,28] is 



q dE 



(2) 



In the modeling process, the density of state (DOS) and the 
Fermi probability function, f(E), are employed. It is notable that 
electrical property of materials from metal to semiconductor is 
changing by the gradient of DOS^) near the Dirac point [29]. 



DOS(£) = 



1 



3na c _ c t 



E 2 - 



■d(E), 



v 2 j 



(3) 



where «c-C = 1-42 A is the carbon-carbon bond length, 
t = 2.7 eV is the nearest neighbor C-C tight binding overlap 
energy, 



g \N + \ 3 



p: integer number , 



is the energy band gap of graphene monolayer, and the Fermi 
probability function f(E) is defined as [30] 
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The integral of these two values with respect to E gives the 
carrier concentration equation as 



n = jlK>S(E)f(E)dE. 



(5) 



By replacing the DOS and f(E) in Equation 3, the carrier 
concentration in the non-parabolic region is defined as 



n = 

Jo 
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exp 
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d(E) (6) 



dQ 
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(i) 



Now the quantum capacitance can be calculated as: 



in which dQ = e-dn is the charge measured in coulombs, e is the 
electron charge, and n is the intrinsic carrier concentration of 
graphene. By substitution of the applied voltage 8V= dEle to 
the device we obtain 



3naQ_Qt 



exp 



E-Ep 



+ 1 (7) 
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The equation provides a quantitative description of the graphene 
quantum capacitance in terms of the Fermi velocity [31], carrier 
density, temperature and fundamental physical quantities. 
According to the relationship between energy band structure 
and the graphene potential, the quantum capacitance-voltage 
characteristic of the proposed model is depicted in Figure 3. In 
the figure, the experimental data is also plotted to have a fair 
scale for validating the proposed model [32]. 



10 
9.5 



* Experimental data 
Proposed model 




V n (V) 



Figure 3: The proposed model of quantum capacitance of EGFETs 
based single-layer graphene. 



To get a greater insight into the quantum capacitance of 
graphene -based EGFET devices, a number of important charac- 
teristics of the C-V curve are discussed. To begin with, at the 
Dirac point, the quantum capacitance has a minimum value 
which is close to zero. The other evident aspect in the quantum 
capacitance model is the linear rise of the capacitance with the 
voltage, which is symmetric with respect to the Dirac point. 
Despite of all mentioned vintages of the proposed model, the 
characteristics of the proposed model diverge considerably from 
the experimental data. To ease this error, a new function of E is 
multiplied by the previous model in Equation 7. This function is 
a square root function, which must be symmetric to its own 
origin. Otherwise, it will disrupt the isochronism of the 
proposed model. Equation 8 shows the general form of the 
suggested function multiplied by the presented model in Equa- 
tion 7. The quantum capacitance, C q , is now 



C=e z 



1 



3naQ_Qt 



r E^ 2 



v 2 y 



exp 



+1 



(ctxj£| + |p|), 



(8) 



where a and p are unknown parameters, which need to be 
adjusted properly. Finding the best fitted values for a and P 
requires an optimization technique with an accurate and reli- 
able performance. To this end, ant colony optimization (ACO) 
is used as one of the well-known and efficient metaheuristic 
swarm intelligence-based optimization algorithms. The ACO 
algorithm has several advantages over conventional mathe- 
matic algorithms [33,34]. It is a fast converging algorithm with 
the capability of escaping from local optima in the search space. 
The random values used in the movements of the particles help 
the algorithm to stochastically improve the obtained solution 
during each iteration. 

Ant colony optimization overview 

Ant colony optimization (ACO), which is inspired by the 
foraging behavior of real ant colonies exploring for foods, was 
proposed by Dorigo in 1992 [33,35,36]. ACO mainly imitates 
the team work of an ant colony in finding a food source. If an 
ant finds a food source, it will carry a portion of the food to the 
nest, after performing some evaluations about the size of the 
source. On the way back to the nest, it releases some phero- 
mones, which are known as pheromone trail. The rest of the 
ants in the nest can reach to the food source by tracing the 
remaining pheromones. The same behavior is shown by the rest 
of the ants in returning to the nest from the source. The amount 
of the pheromones deposited on the way is quite dependant on 
the quantity and quality of the food source [37]. The pheromone 
(x t ) is a vaporizable substance and its amount decreases over 
time. Therefore, the path with the highest amount of pheromone 
is the one, which was chosen by the more ants than the other 
paths. The shortness of the path is a priority to the pheromone 
trail and the ants try to find the shortest possible path. In the 
ACO technique, the pheromone trail that represents a better 
solution, is updated consequently and there exist several ant 
colony models in the literature [37-39]. Equation 9 presents the 
location of the k-th ant in the solution space: 



:x gbest ±dx, where t = 1,2,3,... J 



(9) 



where T is the total number of iterations and t denotes the itera- 
tion number, x gbest is the location of the best objective value 
obtained until iteration t, and dx is a random vector generated 
from [-a, a] to determine the allowed variation the ant can have 
from the x gbest with the same dimensions. The length of this 
jump (variation) is obtained from Equation 10 at the end of Vr 
iterations. 



a t+l =0.lxa t 



(10) 
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The direction of the variation from x gbest , which is shown by ± 
in Equation 9 is decided based on the following equation: 

x best =x gbest +(x gbest x0.0l) 

f/(x best )-/(x gbest )<0 ( + )sign (") 



otherwise 



(-) sign 



To simulate the evaporation of the pheromone, Equation 12 is 
presented, and Equation 13 shows the increment of the 
pheromone around the best objective value obtained after each 
iteration. 



T, =0.1XT M 



(12) 



-(0.01x/(x^ est )) (13) 



Optimization of the proposed model 

The aim of the optimization is to find the best values for a and P 
in Equation 8. Therefore, the search space of the problem is a 
2D space, which returns values for a and P at each iteration. 
The number of ants hired as agents in the ACO algorithm is set 
to be 100, which requires a matrix of 100 x 2 to store the a and 
P values for all the ants, at each iteration. To evaluate the solu- 
tions proposed at each iteration, a fitness function is defined as: 



Fitness function = cp z - = ^ C l q (&) - C q (&) 



k=l 



(14) 



where C l ^(k) represents the modelled quantum capacitance 
waveform for particle i, C q (£) is the experimental data of the 
quantum capacitance, and cp z - is the fitness value for the z'-th ant 
in the colony. The chosen fitness function calculates the squared 
error between the proposed model and experimental data; 
hence, the lowest value for the fitness function indicates the 
best solution to the a and p values. The strategy of the ACO 
algorithm for the optimization is shown in Figure 4 as a flow 
chart. The best values obtained for the parameters a and p after 
the optimization process are shown in Table 1 . The fitness value 
that is the best of 30 runs of the algorithm, and the respective 
values for the desired parameters are tabulated. Figure 5 shows 



Start 



assignment of initial location 
to each ant in the colony 



Calculate the fitness* value 
for all the ants 



Update the pheromone 
trail 



Set Iteration = 0 



No 



Iteration = Iteration +1 




use the pheromone trail to return 
objective values for each ant. 



simulate the evaporation and 
strengthening of the pheromone 



Update the pheromone trail 
based on the solutions created 



Figure 4: A flowchart of ACO-based algorithm for optimizing the 
quantum capacitance model. 
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Figure 5: Comparison between the proposed single-layer graphene 


quantum capacitance model, the optimized proposed model and the 


experimental extracted data. 









Table 1 : The best values of the optimized parameters over the 30 runs. 


number of runs maximum iteration number best fitness value 


optimized value for a 


optimized value for (3 


30 10,000 3.289-10" 6 


1.0753 


0.724 
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Table 2: The MAPE value of the optimized proposed single layer graphene quantum capacitance model. 



capacitance vs voltage characteristic 


MAPE value (%) 


accuracy based on MAPE (%) 


optimized proposed model 


2.54 


97.46 


proposed model 


11.82 


88.18 



the proposed single-layer graphene quantum capacitance model, 
the optimized proposed model and the experimental extracted 
data, as graphs. 

To evaluate the quality of the optimized model compared with 
the experimental data, the mean absolute percentage error 
(MAPE) is used as an error evaluation parameter, as shown in 
Equation 15. 



1 n 

MAPE = -V 

n k=l 



CJk)-C q (k) 



CJk) 



(15) 



The results for the MAPE for the proposed model and the opti- 
mized proposed model are shown in Table 2. The accuracy 
based on the obtained MAPE value is also reported as the result 
of subtracting the MAPE from 100 percent. Based on the results 
tabulated, the accuracy of the optimized model is more than 
97%, which is in an acceptable range of accuracy. 

The logarithmic convergence profile of the best fitness value 
obtained is plotted in Figure 6. The graph indicates that the 
algorithm converges to the optimized values with an acceptable 
convergence speed after around 1500 iterations. 

It is apparent that there is a favorable agreement between the 
optimized proposed model of graphene-based EGFETs device 




800 1000 1200 1400 

Generation 

Figure 6: The convergence profile of the optimization of the proposed 
model using ACO technique. 



and experimental result. It can be concluded that, the presented 
model can be applied as a powerful tool to optimize the 
graphene-based EGFETs device performance. 

Conclusion 

Graphene as a 2D sheet of hexagonally arranged carbon atoms 
exhibits amazing carrier transport properties and a high sensi- 
tivity at the single-molecule level, which makes it a promising 
material for nanoscale devices. According to the graphene 
structure, it can satisfy the major requirements of a channel in 
electrolyte-gated transistor (EGFET) devices due to its ballistic 
transport, high conductivity, and strong mechanical and elas- 
ticity properties. An analytical modeling of the graphene capac- 
itance as a major characteristic of EGFET is studied in this 
paper and the electrical circuit of the device is discussed. An 
EGFET based structure is employed as a platform and the 
graphene capacitance is studied. In order to enhance the accu- 
racy of the proposed model, an ant colony optimization (ACO) 
algorithm is implemented and we obtained acceptable results 
with more than 97% of accuracy. Finally, for the purpose of 
verification, the C—V characteristic of the optimized model is 
investigated with an existing experimental study and shows an 
acceptable agreement. This paper demonstrates how the opti- 
mized model can be used to predict the capacitance variation of 
graphene in graphene-based devices. 

Acknowledgements 

The authors would like to acknowledge the financial support 
from Research University grant of the Ministry of Education 
Malaysia under project number Q.J1 30000.2523. 04H99. Also 
thanks to the Research Management Centre (RMC) of Univer- 
siti Teknologi Malaysia (UTM) for providing an excellent 
research environment in which this work was completed. 

References 

1. Gibertini, M.; Singha, A.; Pellegrini, V.; Polini, M.; Vignale, G.; 
Pinczuk, A.; Pfeiffer, L. N.; West, K. W. Phys. Rev. B 2009, 79, 241406. 
doi: 1 0. 1 1 03/PhysRevB.79.241 406 

2. Zarbo, L. P.; Nikolic, B. K. EPL 2007, 80, 47001. 
doi:10.1209/0295-5075/80/47001 

3. Burgess, C. P.; Dolan, B. P. Phys. Rev. B 2007, 76, 1 13406. 
doi: 1 0. 1 1 03/PhysRevB.76. 1 1 3406 

4. Wang, X.; Xie, W.; Xu, J.-B. Adv. Mater., in press, 
doi: 1 0. 1 002/adma.201 306041 



608 



Beilstein J. Nanotechnol. 2014, 5, 603-609. 



5. Zou, Y.; Li, F.; Zhu, Z. H.; Zhao, M. W.; Xu, X. G.; Su, X. Y. 

Eur. Phys. J. 82011, 81, 475-479. doi:1 0.1 140/epjb/e201 1-20225-8 

6. Chai, Y.; Gong, J.; Zhang, K.; Chan, P. C. H.; Yuen, M. M. F. 
Nanotechnology 2007, 18, 355709. 

doi : 1 0. 1 088/0957-4484/1 8/35/355709 

7. Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; 
Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Science 2004, 306, 
666-669. doi : 1 0. 1 1 26/science. 1 1 02896 

8. Dong, X.; Zhao, X.; Wang, L; Huang, W. Curr. Phys. Chem. 2013, 3, 
29 1 -30 1 . doi : 1 0.2 1 74/1 8779468 1 1 303030006 

9. Ragheb, T.; Massoud, Y. On the modeling of resistance in graphene 
nanoribbon (GNR) for future interconnect applications. In IEEE/ACM 
International Conference on Computer-Aided Design, 2008. ICCAD 
2008, San Jose, CA, Nov 10-13, 2008; IEEE, 2008; pp 593-597. 
doi:10.1 109/ICCAD.2008.4681637 

10. Adam, S.; Das Sarma, S. Phys. Rev. B 2008, 77, 115436. 
doi: 1 0. 1 1 03/PhysRevB.77. 1 1 5436 

11. Hill, E. W.; Geim, A. K.; Novoselov, K.; Schedin, F.; Blake, P. 
IEEE Trans. Magn. 2006, 42, 2694-2696. 

doi:10.1 109/TMAG.2006.878852 

12. Grassi, R.; Gnudi, A.; Gnani, E.; Reggiani, S.; Baccarani, G. Graphene 
Nanoribbons FETs for High-Performance Logic Applications: 
Perspectives and Challenges. In 9th International Conference on 
Solid-State and Integrated-Circuit Technology, 2008. ICSICT 2008, 
Solid-State and Integrated-Circuit Technology, Beijing, China, Oct 
20-23, 2008; IEEE, 2008; pp 365-368. 

doi:10.1 109/ICSICT.2008.4734555 

13. Liang, G.; Neophytou, N.; Nikonov, D. E.; Lundstrom, M. S. 
IEEE Trans. Electron Devices 2007, 54, 677-682. 
doi:10.1109/TED.2007.891872 

14. Bolotin, K. I.; Sikes, K. J.; Jiang, Z.; Klima, M.; Fudenberg, G.; Hone, J.; 
Kim, P.; Stormer, H. L. Solid State Commun. 2008, 146, 351-355. 
doi:10.1016/j.ssc.2008.02.024 

15. Chen, D.; Tang, L.; Li, J. Chem. Soc. Rev. 2010, 39, 3157-3180. 
doi:10.1039/b923596e 

16. Novoselov, K. S.; Fal'ko, V. I.; Colombo, L.; Gellert, P. R.; 
Schwab, M. G.; Kim, K. Nature 2012, 490, 192-200. 
doi:10.1038/nature11458 

17. Neto, A. C; Geim, A. New Sci. 2012, 2012, 214. 

18. Han, M. Y.; Ozyilmaz, B.; Zhang, Y.; Kim, P. Phys. Rev. Lett. 2007, 98, 
206805. doi:10.1103/PhysRevLett.98.206805 

19. Yang, L; Park, C.-H.; Son, Y.-W.; Cohen, M. L; Louie, S. G. 
Phys. Rev. Lett. 2007, 99, 186801. 

doi: 1 0. 1 1 03/PhysRevLett.99. 1 86801 

20. Novoselov, K. S.; Jiang, D.; Schedin, F.; Booth, T. J.; Khotkevich, V. V.; 
Morozov, S. V.; Geim, A. K. Proc. Natl. Acad. Sci. U. S. A. 2005, 102, 

1 0451-1 0453. doi: 1 0. 1 073/pnas.05028481 02 

21. Ouyang, Y.; Yoon, Y.; Guo, J. IEEE Trans. Electron Devices 2007, 54, 
2223-2231. doi:10.1 109/TED.2007.902692 

22. Heller, I.; Chatoor, S.; Mannik, J.; Zevenbergen, M. A. G.; Dekker, C; 
Lemay, S. G. J. Am. Chem. Soc. 2010, 132, 17149-17156. 
doi:10.1021/ja104850n 

23. Ye, J.; Craciun, M. F.; Koshino, M.; Russo, S.; Inoue, S.; Yuan, H.; 
Shimotani, H.; Morpurgo, A. F.; Iwasa, Y. Proc. Natl. Acad. Sci. U. S. A. 
2011, 108, 13002-13006. doi:10.1073/pnas.1018388108 

24. Warner, J. H.; Rummeli, M. H.; Ge, L.; Gemming, T.; Montanari, B.; 
Harrison, N. M.; Buchner, B.; Briggs, G. A. D. Nat. Nanotechnol. 2009, 
4, 500-504. doi:10.1038/nnano.2009.194 



25.Sorgenfrei, S.; Chiu, C.-y.; Gonzalez, R. L., Jr.; Yu, Y.-J.; Kim, P.; 
Nuckolls, C; Shepard, K. L. Nat. Nanotechnol. 2011, 6, 126-132. 
doi: 1 0. 1 038/nnano.201 0.275 

26. O'Connor, I.; Junchen, L.; Gaffiot, F.; Pregaldiny, F.; Lallement, C; 
Maneux, C; Goguet, J.; Fregonese, S.; Zimmer, T.; Anghel, L.; 
Trong-Trinh, D.; Leveugle, R. IEEE Trans. Circuits Syst, /2007, 54, 
2365-2379. doi:10.1 109/TCSI.2007.907835 

27. John, D. L.; Castro, L. C; Pulfrey, D. L. J. Appl. Phys. 2004, 96, 
5180-5184. doi:10.1063/1. 1803614 

28. Chu, L; Xue, Q.; Zhang, T.; Ling, C. J. Phys. Chem. C2011, 115, 
1 521 7-1 5224. doi:1 0.1 021/jp2030768 

29.Schaper, A. K.; Wang, M. S.; Xu, Z.; Bando, Y.; Golberg, D. Nano Lett. 
2011, 11, 3295-3300. doi:10.1021/nl201655c 

30. Karimi, F. A. H.; Ahmadi, M. T.; Rahmani, M.; Akbari, E.; Kiani, M. J.; 
Khalid, M. Sci. Adv. Mater. 2012, 4, 1142-1147. 
doi:10.1166/sam.2012.1405 

31. Venkatesan, B. M.; Bashir, R. Nat. Nanotechnol. 2011, 6, 615-624. 
doi:10.1038/nnano.2011.129 

32. Xia, J.; Chen, F.; Li, J.; Tao, N. Nat. Nanotechnol. 2009, 4, 505-509. 
doi:10.1038/nnano.2009.177 

33. Dorigo, M.; Birattari, M.; Stutzle, T. IEEE Comput. Intell. Mag. 2006, 1, 
28-39. doi:10.1109/MCI.2006.329691 

34. Dorigo, M.; Blum, C. Theor. Comp. Sci. 2005, 344, 243-278. 
doi:10.1016/j.tcs.2005.05.020 

35. Dorigo, M.; Gambardella, L. M. IEEE Trans. Evol. Comput. 1997, 1, 
53-66. doi:10.1 109/4235.585892 

36.Shen, Q.; Jiang, J.-H.; Tao, J.-c; Shen, G.-l.; Yu, R.-Q. 

J. Chem. Inf. Model. 2005, 45, 1024-1029. doi:10.1021/ci049610z 
37.Socha, K.; Dorigo, M. Eur. J. Oper. Res. 2008, 185, 1155-1173. 

doi: 1 0. 1 01 6/j.ejor.2006.06.046 

38. Dorigo, M.; Maniezzo, V.; Colorni, A. 

IEEE Trans. Syst. Man Cybern., Part B 1996, 26, 29-41. 
doi:10.1 109/3477.484436 

39. Rahmani, R.; Yusof, R.; Seyedmahmoudian, M.; Mekhilef, S. 
J. Wind End. Ind. Aerodyn. 2013, 123, Part A, 163-170. 
doi: 1 0. 1 01 6/j.jweia.201 3. 1 0.004 



License and Terms 

This is an Open Access article under the terms of the 
Creative Commons Attribution License 
( http ://creativecommons . org/licenses/by/2 . 0 ), which 
permits unrestricted use, distribution, and reproduction in 
any medium, provided the original work is properly cited. 

The license is subject to the Beilstein Journal of 
Nanotechnology terms and conditions: 
( http://www.beilstein-journals.org/bjnano ) 

The definitive version of this article is the electronic one 

which can be found at: 

doi:10.3762/bjnano.5.71 



609 



